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Abstract. - We measure the fluctuations of the position of a Brownian particle confined by 
an optical trap in an aging gelatin droplet after a fast quench. Its linear response to an external 
perturbation is also measured. We compute the spontaneous heat flux from the particle to the bath 
due to the nonequilibrium formation of the gel. We show that the mean heat flux is quantitatively 
related to the violation of the equilibrium fluctuation-dissipation theorem as a measure of the 
broken detailed balance during the aging process. 
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\ Introduction. — The study of the statistical prop- 
erties of fluctuations in out-of-equilibrium systems is a 
topic of great current interest. For instance, it is im- 
portant at micro- and nano- scales where fluctuations 
and off-equilibrium conditions are common. In particu- 
lar, the fluctuation theorem (FT) [T] and the generaliza- 
tion of the fluctuation-dissipation theorem (FDT) [5H5] 
to processes away from thermal equilibrium, are impor- 
tant theoretical results relevant to describe energy ex- 
changes. A large number of studies has been devoted to 
non-equilibrium steady states and transient states in con- 
tact with a single heat bath. Nowadays, these systems 
begin to be rather well understood both theoretically and 
experimentally (see ref. [7J for a short review). On the 
other hand, non-stationary states slowly relaxing towards 
thermal equilibrium are important for applications. For 
example, glassy systems age in time driven by spontaneous 
heat fluxes to the environment, after being prepared in a 
mctastable configuration [8]. Although several extensions 
of the FT [SHE] and generalized fluctuation-dissipation re- 
lations (GFDRs) PH3S1 have been formulated for relaxing 
systems, the comparison between theory and experiments 
still lacks a clear interpretation. 

In the present letter we study experimentally the fluc- 
tuations, the linear response and the heat flux of a non- 
stationary system: an aging gel during the sol-gel transi- 
tion. The rheological properties and the fluctuations dur- 
ing this aging dynamics can be characterized by a single 
relevant degree of freedom of a local probe. Specifically, 
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we focus, within the context of GFDRs, on the dynam- 
ics of a micron-sized particle embedded in the aging gel, 
which acts as a non-equilibrium bath. The particle is the 
probe, which measures the out-of-equilibrium properties 
of the gel. Our aim is to show in a transparent way the 
link between the violation of the equilibrium FDT and 
the total entropy production in an experiment performed 
in an aging system. 

Description of the system. — In the experiment, a 
silica bead of radius r = 1 /im is used as a probe inside a 
thermoreversible gel (gelatin) slowly relaxing towards its 
solid- like state (gel), after a very fast quench, from above 
to below the gelation temperature T ge i. Above T ge i an 
aqueous gelatin solution is in a viscous liquid phase (sol), 
whereas below T ge i the formation of a network of cross- 
linked triple helices leads to an elastic solid-like phase (gel) 
|13j . This relaxational out-of-equilibrium regime is probed 
by measuring the Brownian fluctuations and the linear 
response function of the bead position. 

Specifically, the gel used in the experiment is an aque- 
ous gelatin solution (type-B pig skin) at a concentration 
of 10 wt %, prepared following the usual protocol [T3JU5]. 
For this sample T ge i — 29° C. The solution fills a transpar- 
ent cell kept at constant temperature (To = 26 ± 0.05°C 
< T ge i) by a Peltier element. See the sketch of the ex- 
perimental setup on the left panel of fig. da). Thus, the 
solution inside the cell is in the solid-like phase. A sil- 
ica bead, of radius r = 1 /jm, embedded the gelatin so- 
lution is placed in the focal position of a tightly focused 
laser beam (A = 980 nm) at a power of 20 mW. At this 
power the laser produces on the particle an elastic force 
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Fig. 1: (a) Sketch of the experimental setup used to melt a 
micron-sized volume of gel by locally heating the bulk at To = 
38° C with a focused laser beam, (b) Diagram of the quench of 
the liquid droplet from T = 38° C to T = 27.2° C in less than 
1 ms. Lower panel: position fluctuations of a particle optically 
trapped inside the liquid before and after the quench. 



of stiffness k = 2.9 pN//zm. Because of light absorption, 
the temperature of the trapped particle is T = 27.2°C, 
which is still smaller than T ge i. Hence, the bead is inside 
the solid gel in the beam focus at a distance ft, = 25 /Ltm 
from the cell wall. Starting from this condition, the laser 
power is increased to 200 mW and the local temperature 
around the focus rises to 38° C > T ge i [16]. As a result 
the gel melts and a liquid droplet of radius a = 5 /jm, 
is formed around the trapped bead inside the the solid 
gel bulk, as depicted on the right panel of fig. QJa) . Af- 
ter 180 s, the laser power is suddenly decreased again to 
20 mW and the temperature is homogenized by heat dif- 
fusion into the gel bulk. Since the thermal diffusivity of 
water is k = 1.4 x 10~ 7 m 2 s _1 , the heat diffuses in a time 
a 2 / k ~ 0.2 ms and the droplet is very efficiently quenched 
to the final temperature T = 27.2°C < T gei . This tem- 
perature remains constant after the quench due the Peltier 
element at To = 26 ± 0.05°C. In addition, we check that 
the trap stiffness k — 2.9 pN//xm remains constant as well 
throughout the gelation process. At T the liquid inside the 
droplet completely solidifies in about 1 hour and the parti- 
cle, trapped in the center of the drop by the focused beam, 
is a probe of this relaxation dynamics. The quenches are 
repeated up to 60 times in order to compute the proper 
ensemble averages of the measured quantities |17j . 

Active microrheology. — Immediately after the 
quench we record the time evolution of the x position [sec 



fig. []Jb)] °f the trapped particle measured by a position 
sensitive detector whose output is sampled at 8 kHz and 
acquired by a computer. The resolution of the measure- 
ment of x is better than 1 nm [I51IT5] . 

Due to the aging process of the gelatin droplet, the time 
series of x shown in fig. E^b) is non-stationary for time 
t > after the quench. In order to characterize this aging 
time evolution in the linear response regime, we perform 
active microrheology. For this purpose we apply an ex- 
ternal driving to the trapped bead. This is realized by 
displacing the position xq of the optical trap, using a pre- 
cise actuator, which deflects the trapping laser beam. A 
sinusoidal modulation of xq results in a time-dependent 
force F = kxQ on the bead: F(f,t) — Fo sin(27r / t) is 
applied at different driving frequencies 0.2 Hz < / < 5 Hz 
and fixed amplitude Fq = 87 fN. By measuring in the 
time interval [t, t + At] the response R(f, t) of the par- 
ticle position to F(f,t), we obtain the shear modulus 
G(f,t) = G'(f,t) - iG"(f,t) of the droplet at different 
times t after the quench, where G' is the storage and G" 
the loss modulus [T5] (see also eqJSJ. 

Because of the relatively fast aging process due to the 
smallness of the droplet, the spectral analysis involved in 
the calculation of G' and G" must be carried out over a 
short time window [i, t + At] for each aging time t. In 
the following we set At = 15 tQ. As the response R(f,t) 
converges more rapidly than fluctuations, the values of G 
are averaged only on 20 quenches for each value of f. In 
figs. [5Ja) and [2fb) we plot the time evolution of G' and 
G" measured at / = 5 Hz during the first 20 minutes after 
the quench. We can identify three different regimes: 

I. For s < t < 60 s, G' is completely negligible 
whereas G" increases continuously in time by a fac- 
tor of almost two. In this aging regime, hightlightcd 
in fig. |5fb), the gelatin droplet behaves as a purely 
viscous liquid even when the final temperature of the 
quench is below T ge i. 

II. For 60 s < t < 200 s, G" continuously increases and 
G' starts to increase slowly, as plotted in fig. HJb). 
This result shows that the droplet is in a transient 
regime towards the sol-gel transition where the gel 
network is not completely formed. Indeed, G' is still 
much smaller and it increases slower than G" . 

III. For 200 s < t < 1200 s, both G' and G" reach a 
logarithmic growth ~ logi where G' increases faster 
than G". This growth as t increases is similar to 
that reported in macroscopic bulk measurements |13l 
[14l[20] but taking place much faster because of the 
smallness of the gelatin droplet. Then the system is 
actually undergoing gelation providing evidence that 
the percolating gel network is already formed. 



lr This value is large enough to resolve the frequencies of the ap- 
plied driving and at the same time it is short enough to avoid a 
pronounced time evolution of the nonstationary signal x. 
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Fig. 2: (a) Time evolution of G' and G" of the gelatin droplet 
after the quench, measured at / = 5 Hz. (b) Expanded view 
during the first 120 s. The error bars represent the standard 
deviation computed from 20 independent quenches. 
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Fig. 3: (a) Storage and (b) loss modulus mea- 
sured at different times after the quench: t = 
60 (*), 120 (o), 180 (□), 300 (O), 600 (<), 900 (>) and 1200 s (A). 



In figs. [3]Ja) and l^b) we plot the low-frequency spec- 
trum (0.2 Hz < / < 5.0 Hz) of G' and G" , respectively, 
at different times after the quench. In figs. the val- 
ues of G' in regime I (0 s < t < 60 s) are not reported 
as their mean values are close to zero and smaller than 
the statistical errors. Besides, in regime I G" is approxi- 
mately proportional to / [dashed line in fig.[3jb)], and the 
dynamic viscosity of the droplet G" /(2nf) is frequency- 
independent. Therefore, at this time scales we check that 
the fluid inside the droplet actually behaves as a Newto- 
nian fluid with a vanishing relaxation time. In contrast, 
for t larger than 60 s, G' begins to grow. In fig. [3ja) we 
plot as dashed lines the curves oc / and oc f 2 . As / — > 0, 
G' displays an intermediate behavior between these two 
reference curves. Then, we conclude that G'(f,t) — > as 
/ ->■ in regime II (60 s < t < 200 s). On the other hand, 
in fig. [3jb) we observe that G" is roughly proportional to 
/. This behavior of G can be approximately described by 
the Maxwell model for a viscoelastic fluid with a single 
relevant relaxation time tq and zero-shear viscosity rjo 



G(f,t) 



l + 2mfT (t)' 



(1) 



Based on the Maxwell model, we estimate To(t) and the 
zero-shear viscous drag coefficient 70 (t) = 6Trrrj (t) at ag- 
ing time t (see ref. [T5] for details) 

As expected, To is very small (« 2 ms) during the first 
60 s after the quench. Therefore, in regime I, the sol 
droplet can be regarded as a Newtonian fluid. Instead, for 
60 s < t < 200 s the mean value of To becomes four times 
the value found for the regime I. The values of To must be 
compared with the viscous relaxation time of the particle 
inside the optical trap: t^ = 70/fc. Since tq/tu < 0.05 
even in regime II, the viscoelastic memory effects of the 
droplet taking place during To do not significantly affect 
the Brownian motion of the particle during the first 200 s. 
Finally, for regime III (200 s < t < 1200 s) G" and mainly 
G' exhibit a complex frequency dependence. Unlike the 
Maxwellian behavior [G 1 — > as / — > 0), in this regime 
the value of G' seems to remain non-zero as / —> 0, as 



shown in fig. [3Ja) . This result is in agreement with bulk 
measurements |13II14] that report a non-zero storage mod- 
ulus of constant value at very low frequencies. This low- 
frequency behavior corresponds to that of an elastic solid 
with multiple relaxation time-scales. 

Heat flux. — We now focus on the spontaneous energy 
fluctuations (F = 0) of the Brownian particle inside the 
aging droplet for the very first 200 s after the quench. 
Since the motion of the trapped particle is overdamped, 
its total energy at time t is 



U(t) 



1 



kx 2 + U s tored{t), 



(2) 



where U st0 red(t) is the energy stored by the surrounding 
gelatin chains in the droplet until time t. By considering 
a non-Markovian overdamped Langevin equation for the 
dynamics of x and eq. (HJ as a model for the corresponding 
memory kernel [19] . the mean energy stored by the bath 
at time t can be estimated in this time interval (aging 
regimes I and II) by [T7] 



(U s tored(t)} — 



(3) 



As 2lM < 0.05 for the first t < 200s, we deduce from 
eq. ([3]) that the energy stored by the bath constitutes in 
average less than 5% of the total potential energy of the 
particle (see also ref. [H]). Hence, the instantaneous value 
of the total potential energy at time t is given with an error 
smaller than 5% by Ut = kx 2 /2. As there is no external 
force acting on the particle, the heat exchanged between 
the particle and the bath during the time interval [t, t + r] 
is equal to the variation AJ7{ jT = U t + T — U t of the total 
energy of the particl^l (see ref. [H]). Specifically: Q t , T = 
AUt, T = k(x 2 +T —xf)/2, where a positive (negative) value 
of Qt.r represents a heat fluctuation from (to) the bath to 



2 This is not valid during the aging regime III because in such a 
case we have to take into account the particle history between t and 
t + r determined by the exact form of the shear modulus G(f,t). 
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(from) the bead. Then, the mean heat transferred during 
[t, t + t] can be written as 



(Qt,r) = ^[^(i + r) 2 



(4) 



where cr x (t) 2 is the ensemble variance of the spontaneous 
fluctuations of x at time t. In ref. [T5] we show that the 
fluctuations of x are Gaussian and <J 2 (t) is a monotoni- 
cally decreasing function of t. This time evolution is plot- 
ted in fig. [Ha) for the three different aging regimes iden- 
tified by active microrheology. During the first « 20 s, 
Ca;(i) 2 > ksT/k, which can be interpreted as the result 
of both thermal and non-thermal stochastic forces on the 
particle due to the transient formation of the gel network, 
as schematized in the inset of fig.[Ha). The dimensionlcss 
heat q t . T = Qt,r 1 1 {ksT), computed using eq. (|4|), is plotted 
in fig.[DJb). Because of the relaxation of a x (t) 2 , we observe 
the existence of a mean heat flux from the particle to the 
surroundings on the time scale r. The absolute value of 
the mean heat increases as the measurement time r in- 
creases and the maximum value \(qt, T )\ ~ 1 takes place 
at t = s. Note that at thermal equilibrium the mean 
heat would be (<?t, T ) = for all t and t. Non-negligible 
values of the mean heat persist for several seconds after 
the quench. Nevertheless, as t increases, |(<7t. T )| decreases 
becoming experimentally undetectable for i > 20 s. The 
non-vanishing mean heat flux (qt, T ) is a signature of the 
broken detailed balance of the particle dynamics due to 
the assembling gelatin chains. However, as the bath is un- 
dergoing aging the time scales of the dynamics slow down. 
Then the rate at which the heat flows from the system 
to the environement becomes undetectable by the Brow- 
nian probe, which results in an apparent equilibrium-like 
behavior for t > 20 s. In order to check for possible experi- 
mental artifacts, we repeated the same quench experiment 
using a Newtonian fluid (glycerol 60 wt % in water). At 
this concentration the glycerol solution has the same vis- 
cosity of the initial sol phase of gelatin (9 x 10~ 3 Pa s) but 
with no sol-gel transition in the same temperature range. 
For this fluid, the measured mean heat flux is plotted for 
comparison in fig. 3J In this case we verify that the parti- 
cle quickly equilibrates with the bath after the quench so 
that (qt,r) = within the experimental accuracy. 

In ref. [T5] we show that the heat fluctuations verify the 
FT. In analogy with a system kept in a non-equilibrium 
steady state in contact with two reservoirs at unequal tem- 
peratures [2"2"H2"9"j . the quantity 
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<J x (t) 2 <r x (t + r) s 



Qi 



(5) 



represents the total entropy produced by the breakdown 
of the time-reversal symmetry due to the nonstationarity 
of the bath after the quench. By taking the ensemble av- 
erage of eq. ([5]) and using eq. (U]), one finds that the mean 
heat (Qt,r) actually quantifies the mean total entropy pro- 
duction between time t and t+r: (ASt, T ) oc 
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Fig. 4: (a) Time evolution of the variance of the particle 
position after the quench. The insets depict the possible con- 
figurations of the bead during the sol-gel transition of the sur- 
rounding gelatin, (b) Mean value of the normalized heat at 
different times t after the quench in gelatin, for different values 
of the time lag: r = 0.25 s (o), 2.5 s (□), 10 s (O), 30 s (v). 
The dashed line correspond to the normalized mean heat for 
the quench in glycerol at r = 30 s. Inset: Comparison between 
the right-hand side (RHS) and left-hand side (LHS) of eq. |TT}, 
normalized by 2/csT/fc, computed over a time window [t, t+At] 
for At = 15 s. 



Generalized fluctuation-dissipation relation. 

We now discuss the relation between entropy produc- 
tion, spontaneous fluctuations and linear response of the 
trapped particle in this out-of-equilibrium aging bath. It 
should be expected that there is a strong connection be- 
tween the deviations from the equilibrium FDT and the 
extent of the broken detailed balance, quantified by (Qt.r)- 
If the entropy production takes place at a very slow rate, 
then the particle in the aging bath should exhibit an 
equilibrium-like behavior, see the experiments on FDT in 
colloidal glasses of refs. [TSIED] ■ On the other hand, if 
the rate is fast enough, then currents are non-negligible 
and such deviations should be significant like in the non- 
equilibrium steady state experiments of refs. - In 
the present case of gelatin, we found that \{Qt, T ) \ is com- 
parable to fcgT during the first 20 s after the quench. 
The rate of this heat flux slows down as the gelatin ages. 
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Fig. 5: Passive and active spectra computed on the interval 
[t, t + At] at different times t after the quench: (a) [0 s, 15 s], 
the dotted circle points out the low-frequency deviation from 
eq. ([Sj) ; (b) [30 s, 45 s], the arrow indicates the position of the 
frequency f c ; (c) [75 s, 90s]; and (d) [1200 s, 1215 s]. The 
dotted line corresponds the Lorentzian curve obtained using 
eq. © (G' = 0). 

Therefore, a significant deviation of the equilibrium FDT 
is expected to occur for the position of the trapped par- 
ticle during the aging regime I of the gelatin droplet. In 
contrast, in regimes II and III the FDT must be appar- 
ently satisfied with the temperature T of the thermal bath. 
Indeed, we show in the following that our experimental re- 
sults display such a kind of behavior. 

As we performed active microrheology at a fixed ex- 
citation frequency /, it is more convenient to study the 
fluctuations and linear response in the frequency domain. 
The Fourier transform R of the linear response function 
of x to the perturbative time-dependent force F is related 
to the shear modulus G of the gelatin droplet by 

k U' t)= *«rG*tf,t) + k ' (6) 

where G* is the complex conjugate of G. Then, the imag- 
inary part of eq. ([6]) is 



lm{R(f,t)} 



(7) 



[fc + GnrG'if, t)] 2 + [6nrG"(f, t)] 2 

At thermal equilibrium, the quantity given in eq. © docs 
not depend on time and is related to the power spectral 
density (\x(f,t)\ 2 ) of the spontaneous fluctuations of x in 
the absence of an external force (F = 0). This relation is 
provided by the FDT and reads 



(\x(f,t)\ 2 ) = ^h-n{R(f,t)}. 



(8) 



During the non-equilibrium gelation of the gelatin droplet, 
cq(5] does not necessarily hold because as we have already 
seen R is a function of both t and /. Thus we study 
the relation between the passive spectrum (\x(f,t)\ 2 ) and 
the active spectrum 2fcsTIm{i?(/, t)}/(jrf), as a function 
of t. For this purpose we compute independently these 
quantities over a time interval [t, t+At) of length At = 15 s 
at time t after the quench in the different aging regimes. 

In fig.[5ja) we plot (\x(f,t)\ 2 ) for t = s, where a strong 
mean heat flux takes place (regime I). On the other hand, 
taking into account the negligible values of To in this 
regime, Im{i?} can be written as 



Im{R(f,t)} = 



1 



/ 



27r 7 o(i) P + fc{tf 



(9) 



where f c (t) is related to the viscous relaxation time (t) = 
7o(£)/fc of the trapped particle at time t by: 



fc(t) 



2-KT k {t)' 



(10) 



In fig. (5{a) we also plot the active spectrum 
2kBTIm{R}/(irf) computed using eq. ([7]). This spectrum 
has a Lorentzian profile. As expected, an important 
deviation of eq. is observed for frequencies / < 1 Hz. 
Note that for frequencies / > 1 Hz eq. (0 is verified, 
indicating that the deviation is actually due to the slow 
structural assemblage of the ge{f|. The spectral curves 
of fig. [SJa) provide information on the time scales at 
which this phenomenon perturbs the Brownian probe: 
> 1 s. The extent of the deviation, i.e. the difference 
between the area below (\x(f,t)\ 2 ) minus the area below 
2fc B TIm{i?}/(7r/) in fig.^a), is 



<|i(/,*)i 2 > 



2k B T 



lm{R(f,t)} 



df 



2Q 



t,At 



fc 



(11) 

where Q t At can be interpreted as the heat that the system 
dissipates in average during the time interval [t, t+At] [17]. 
This interpretation is justified because, using eqs. (jTTj) . © 
and (UJ) , it can be easily shown [T7] that indeed Q t At = 

~ St It +At (Qt'.f) dt' with f > 30 s, such that at time 
t' + f the systems has reached equilibrium for all times 
after the quench [see figQJa)]. We used the experimental 
data to compute independently the left-hand side (LHS) 
and the right-hand side (RHS) of eq (JTTJ) and we find that 
they are equal within error bars, as shown in the inset of 
fig. lib). 

Thus similarly to the GFDR for non-equilibrium steady 
states derived in [34] , eq. (fTTj) hightlights the role of heat 
currents (or equivalently the entropy production). In the 
present case, the right-hand side of eq. (| 1 1 [) quantifies the 
heat excess that must irreversibly flow from the particle 
to the bath in order for the system to reach equilibrium 



,! Similar low-frequency violations of the FDT are found in non- 
equilibrium active systems, e.g. biological polymer networks 1331 . 
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at temperature T. Thus, the violation of the FDT ([5J 
can be interpreted as a measure of the heat that must 
be dissipated to reach an equilibrium state. This is in 
close analogy to the power dissipated into the medium to 
keep a non-equilibrium steady state in the GFDR of ref. 
[Mj . In fig. EJb) we plot the passive and active spectra 
in the regime I for the t = 30 s. We check that in this 
case the equilibrium-like cq. (JSJ) is satisfied with the final 
equilibrium temperature T = 27.2°C. This results is in 
good agreement with the fact that the mean heat flux 
becomes undetectable by the Brownian particle for t > 
20 s. 

Fig. 0c) displays the passive power spectrum of x for 
t = 75 s, (regime II) where G' starts to increase. On 
the other hand, taking into account that the fluid can be 
characterized by a single relaxation time To < 0.05Tfc, the 
active term can be approximated with an error smaller 
than 5% by eq. © ■ Using this expression with the exper- 
imental values of 70 and f c , in fig. [Sfc) we plot the active 
term 2£^TIm{i?}/(7r/). In this case, the agreement be- 
tween the left- and the right-hand side of eq. (jSJ) is much 
better than in fig. Hence, we verify that the system 

is relaxing to an equilibrium-like behavior. 

Finally, in fig. [5jd) we plot the passive and active spec- 
tra for t = 1200 s. In this case the large viscoelasticity of 
the gelatin droplet gives rise to the highly non-Lorentzian 
shape of both terms. For comparison, we also plot the 
Lorentzian profile (dotted line) that would be obtained 
using equation (|9|) without taking into account the contri- 
bution of the storage modulus G'. In the absence of an 
accurate model for G in regime III, we only plot the active 
term in the frequency range 0.2 Hz< / < 5 Hz at which 
the sinuosoidal external force was applied to the particle. 
We verify that an equilibrium-like fluctuation-dissipation 
relation (|8) holds for the particle position x because of the 
smallness of the heat fluxes as the gelatin droplet ages. 

Conclusion. — We have experimentally measured the 
fluctuations of the position of a trapped Brownian parti- 
cle in a non-stationary bath, i.e. an aging gel after a 
very fast quench. This simple experiment has allowed us 
to understand several concepts formulated originally for 
non-equilibrium steady states that can be naturally ex- 
tended to non-stationary systems slowly relaxing towards 
an equilibrium state. In particular, we provide evidence 
that the extent of the violation of the equilibrium FDT is 
directly related to the entropy that must be produced in 
order for the system to reach equilibrium. Therefore, as 
the system ages, the relation between fluctuations and lin- 
ear response behaves like in equilibrium. This approach to 
fluctuation-dissipation, alternative to the concept of effec- 
tive temperature [5] , can be exploited to study small heat 
fluxes in experiments on more complex aging systems. 

* * * 
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